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Abstract. We investigate the phase diagram of the Heisenberg antiferromagnet on 
the square lattice with two different nearest-neighbor bonds J and J' (J- J' model) at 
zero temperature. The model exhibits a quantum phase transition at a critical value 
J' c > J between a semi-classically ordered Neel and a magnetically disordered quantum 
paramagnetic phase of valence-bond type, which is driven by local singlet formation on 
J' bonds. We study the influence of spin quantum number s on this phase transition 
by means of a variational mean-field approach, the coupled cluster method, and the 
Lanczos exact-diagonalization technique. We present evidence that the critical value 
J' c increases with growing s according to J' c oc s(s + 1). 
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1. Introduction 

The study of quantum antiferromagnets in low- dimensional systems has attracted 
much attention in recent years, both theoretically and experimentally. In particular, 
quantum phase transitions are in the focus of interest, see e.g. Refs. [U El 0]. For 
these zero-temperature transitions thermal fluctuations are irrelevant and the transition 
between different quantum phases (e.g. between magnetically ordered and disordered 
phases) is driven purely by quantum fluctuations. For the quantum spin Heisenberg 
antiferromagnet (HAFM) on two-dimensional lattices the interplay of interactions and 
fluctuations is well balanced and the existence of semi-classical magnetic long-range 
order depends on the degree of competition between bonds jHEj. Competition between 
bonds in spin systems may appear as frustration, which is present in classical as well 
as in quantum spin systems. In quantum systems also a direct competition between 
bonds exists which may lead to local singlet formation on certain antiferromagnetic 
bonds (or plaquettes of four spins) if these bonds are increased in strength. By tuning 
the degree of competition zero-temperature order-disorder phase transitions can be 
realized. The existence of magnetically disordered quantum paramagnetic ground states 
in (quasi-) two- dimensional Heisenberg systems has been recently demonstrated, e.g. for 
SrCu 2 (B0 3 ) 2 0IZ1 and CaV 4 9 jEH- 

A canonical model to study the competition in a frustrated quantum spin HAFM 
is the J\-Ji model on the square lattice, where the frustrating J 2 bonds plus quantum 
fluctuations lead to a second-order transition from Neel ordering to a quantum 
paramagnetic phase, see e.g. [TQl [TTJ [TJJ HS1 E3- A widely studied model describing 
competition without frustration and showing the 'melting' of semi-classical Neel order 
by local singlet formation is the HAFM on the square lattice with two non-equivalent 
nearest-neighbor bonds J and J' (J- J' model) [T01 HZl HH1 HZl H3 I2H 1221 In 
these papers on the J- J' model the extreme quantum case s = 1/2 is considered and the 
competition can be tuned by variation of the exchange bond J'. One finds a second-order 
transition from the quasi-classically Neel ordered phase to a dimerized singlet phase at 
J' c ~ 2.5 . . . 2.9 J. It is argued in Refs. [201 EE! that the quantum phase transition is of 
the same universality class as the thermal phase transition of three-dimensional classical 
Heisenberg model. 

The strength of quantum fluctuations within this model can be varied either by 
anisotropy or by spin quantum number. Indeed its was found in J2Hj for the J- J 1 model 
that the critical J' c for the XY model is significantly larger than for the spin rotationally 
invariant Heisenberg model. The influence of an Ising exchange anisotropy A/ leads also 
to an increase of J' c which is in good approximation proportional to A/ [22! ■ The role 
of the spin quantum number s was not systematically studied. Some results for spin 
models with s = 1 can be found in e.g., [2U 121! 

In the present paper we study the ground state phase transition between a Neel 
ordered phase and a dimerized singlet phase of the J- J' model with spin quantum 
number s = 1/2,1,3/2,2 using a variational mean-field like approach (MFA), the 
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coupled cluster method (CCM) and exact diagonalization (ED) of finite systems. 



2. Model 



We consider the J-J' model on a square lattice, i.e. a HAFM with two kinds 
of antiferromagnetic nearest-neighbor bonds J and J' (see Fig^) described by the 
Hamiltonian 



H 



+ j'Y /Si - Sj , (i) 

<ij>i <ij>2 

where the sums over < ij >i and < ij >2 represent sums over the nearest-neighbor 
bonds, shown in Fig^as dashed and solid lines, respectively. We consider spin operators 

= s(s + 1) of spin quantum number s = 1/2, 1, 3/2, 2. 
Each square-lattice plaquette consists of three J bonds and one J' bond. In what follows 
we set J = 1 and consider J' > 1 as the parameter of the model. In the classical limit 
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Figure 1. Illustration of arrangement of bonds in the J-J' model on the square lattice 
(Eq.Ql): J - dashed lines; J' - solid lines; A and B characterize the two sublattices 
of the classical Neel ground state. 

the ground state is the symmetry breaking Neel state. However, quantum fluctuations 
may lead to a rotationally invariant dimerized valence-bond state for large enough J' . 

3. Variational mean- field like approach 

In this section we use the MFA to calculate the ground-state magnetic order parameter 
of the J-J 1 model (0). For the spin half HAFM this approach has been successfully 
applied to bilayer systems |2Sj, to the isotropic EE] and anisotropic |22j J-J' model 
on the square lattice and on the 1/5 depleted square lattice for CaX^Og pj, but also on 
the strongly frustrated HAFM on the star lattice [26J. In this paper we extent the basic 
ideas of this approach to higher spin quantum numbers s = 1, 3/2, 2. 
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We start with the two uncorrelated mean-field like states, namely the Neel state 
I^MFi) = I + s)\ — s)\ + s)\ — s) . . . and the dimerized rotationally invariant singlet 
product state (valence-bond state) \4>mf 2 ) = n<ij> 2 \{hj}s), where the product runs 
over all J' bonds, cf. ((TJ). \{(i,j)} s ) in \4>mf 2 ) is a singlet state of two spins s, i.e. we have 

IO.W = -M+h)\-h)-\-\)\+h), \{hj} s =i) = ^[i+i)i-i)-|o)|o)+i-i)i+i)], 

\{hj}s=3/2) = ll\ + |)| - |) - I + |)| -|) + |- |)l + \) ~ I " |)l + I)], \{i,j}s=2) = 

jg[\ + 2)| - 2) - | + 1)| - 1) + |0)|0) - | - 1)| + 1) + | - 2)| + 2)]. 

In order to describe the transition between both states we consider for the different 
spin quantum numbers respective uncorrelated product trial states of the form 



I var I 

|^ S=1 ) 
I var I 

|^= 3 / 2 \ 
I var I 



VTT 



l 



l 



+ 2 )l "2 ) - a| -2 )l + 2 ) 



+ 1)| -l)-6i|0)|0) + 62|-1)| + 1) 



c\ + c| 



+ §>l-§>-«*l + 5>l-5> 



I var I 



n 

<«j>2 

n 

<lj>2 

n 

<ij>2 

n - 

^ y/l + d? + d\ + d§ + c£ 

+ 4|0)|0) - d 3 | - 1)1 + 1) + 4| - 2)| + 2) 



(2) 
(3) 



+ c 2 |--)| + -)-c 3 | --)! + -) 



(4) 



+ 2)| -2) -4| + 1)|-1) 



(5) 



where in the two-spin states \n)\m) the first bra vector belongs to site i and the second 
to site j of a J' bond. The trial wave functions depend on the variational parameters 
a; 6i, 62; ci, C2 C3; 4, 4, 4, 4 and interpolate between the valence-bond state \4>mf 2 ) 
realized for a — 1; 61 = 62 = 1; Ci = C2 = C3 = 1; 4=4 = 4 = 4 = 1 and the 
Neel state |0a/Fi) for a = 0; b\ = 6 2 = 0; ci = C2 = C3 = 0; 4 = 4 = 4 = 4 = 0, 
respectively. The ground-state energy per site e s var = (^l ar \H\^l ar ) / N is calculated as 



e var ^ ( a ) 



f_ a + 1(1 + a 2 ) 
~2~ 1 + a 2 



3 (1 



,2\2 



2 4(1 



,2\2 



,s=l 



(6l,62 



26i6 2 + 1+61 



(1-6|) 2 



<ar 3/2 (Cl,C2,C 3 ) 



J' 26i 

~2 l + b\ + 62 2 (1 + b\ + 6 2 ) 2 

J' 3ci + 4cic 2 + 3c 2 c 3 + 1(9 + c\ + c\ + 9c 
2" 



(6) 
(7) 
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e U (4,4,4,4) 



J' 44 + 644 + 644 + 444 + d\ + 4 + 0% + 4^ 2 
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(9) 
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The relevant order parameter describing the Neel order is the sublattice magnetization 
M = <*S«r|s?6Al*S«r>- Usin S Ec l s - ©,©,©,© we obtain 

= (10) 

M^ft.y = - (11) 



1 + 6? + ii 

- 2(i + + tr+^+c1) <i2) 

2 + c/ 2 — rf 2 — 2<i 2 

^*'*-*^ = TT^T^' (13) 




Figure 2. Sublattice magnetization M/s versus J' calculated by the variational 
mean-field like approach (MFA), see text. 

We minimize E = (ty var \H\^ var ) with respect to the variational parameters. As a 
result we obtain an analytic expression for a in the case of s = 1/2, but a set of 2, 3, 
4 coupled nonlinear equations for s = 1, 3/2, 2 to determine the variational parameters. 
As reported in [THJ |2| the sublattice magnetization for s = 1/2 is M S=1 / 2 (J') = 
— (J'/3) 2 for J' < 3 but zero for J' > 3. Furthermore, one can express the ground- 
sate energy as a Landau-type function of M, e»ar' 2 = — §«/' + |(</' — 3)M 2 + |J'M 4 , 
indicating the molecular field-like nature of the approach. For s = 1, 3/2, 2 we have to 
solve the corresponding set of nonlinear equations numerically. We show M(J') in FigHJ 
M(J') vanishes at a critical point J' c = 3 (s = 1/2), j£ = 8 (s = 1), J' c = 15 (s = 3/2), 
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J' c = 24 (s = 2), respectively. The corresponding critical index is the mean-field index 
1/2. 

The sequence of critical points for s = 1/2,..., 2 are precisely described by 
J c (s) = |s(s + l)(z — 1), where z = 4 is the coordination number of the square lattice. 
Although we do not have results for s > 2, we argue that due to the systematic character 
of the MFA approach it seems to be likely that this expression is valid also for s > 2. 

4. Coupled cluster method (CCM) 

We now briefly describe the general CCM formalism, for further details the interested 
reader is referred to Refs. [23 EHl EMI I2E] ■ In order to calculate the many-body ground 
state, we start with a normalized reference or model state |$). We chose the Neel 
state as the reference state |<3>) in order to treat the J- J' model using the CCM. It is 
convenient to perform a rotation of the local axis of the up spins such that all spins in 
the reference state align in the same direction, namely along the negative 2-axis, such 
that we have |<&) = | — s) \ — s)\ — s) \ — s) . . . . We define a set of multi-spin creation 
operators Cf = sf , sfsf , sfsfs^ , . . . . With this definition of the Cf we have 
(<&\Cf = = C7|$), where Cj is the Hermitian adjoint of Cf. The CCM ket and bra 
ground states are then given by 

|*> = e 5 |$>, S = Y,SiCf, (14) 
(|r| = ($|SV 5 , S = 1 + J2^i°i- ( 15 ) 

The correlation operators S and S contain the correlation coefficients Si and Sj which 
have to be determined. Using the Schrodinger equation, H\^) = E\ty), we can now 
write the ground-state energy as E = ($|e" 5 ife s |$) . The sublattice magnetization is 
given by M = -(^|sf |\&) . 

In order to determine the correlation coefficients Si and Si, we to require that the 
expectation value H = {^>\H\^f) is a minimum with respect to Si and Si. In case we 
would be able to take into account all possible multispin configurations in the correlation 
operators S and S the CCM formalism is exact. However, for the considered quantum 
spin model we have to use approximation schemes to truncate the expansion of S and S 
in the Eqs. (114)1 and 1)15)) . As in Refs. |2T ) 128 ) 129). we use the SUBn— n approximation 
scheme, where we include only n spin flips in all configurations (or lattice animals in the 
language of graph theory) which span a range of no more than n adjacent lattice sites. 
Note that this approximation for s = 1/2 is equivalent to the LSUBn approximation 
)27) 123 EH] • Since the approximation becomes exact in the limit n — > oo is useful to 
extrapolate the 'raw' CCM-SUBn— n results to the limit n — ► oo. Although an exact 
scaling theory for SUB results is not known, there is empirical experience |18ll2~Tll2*Hll2Ti] 
how the order parameter for antiferromagnetic spin models scales with n. In accordance 
with those findings we use M{n) = M(oo)+ai(l/n)+a2(l/n) 2 to extrapolate to n — > oo. 
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We note that we take a value of M(oo) tending to zero to indicate the critical point 
J' c (see FigJ3J). The values for J' c obtained by extrapolation of the SUBn— n results for 
M are, however, found to be slightly too large fHJE2j- Therefore, it is more favorable 
to consider the inflection points of the M(J') curve for the SUBn— n approximation, 
assuming that the true M(J') curve will have a negative curvature up to the critical 
point. For increasing n the inflection point J' in j approaches the critical point J' c . The 
inflection points for the SUBn— n approximation again we can extrapolate to the limit 
n — > oo using J' in Jn) = J' in f(oo) + &i(l/n) + &2(l/ n ) 2 and interpret J' in j(oo) as the 
critical value J' c . 




Figure 3. Sublattice magnetization M/s versus J' for spin quantum number s = 1 
using coupled cluster method (CCM), see text. 

In principle it is possible to apply the CCM for arbitrary spin quantum number 
s. However, within the used SUBn— n approximation scheme for higher s additional 
problems appear, namely (i) the number of fundamental configurations (lattice animals) 
increases with s, which makes the calculations on a certain level of approximation n more 
ambitious and (ii) the total number of basis states grows drastically with s according to 
s N and as a consequence the SUBn— n approximation may become less reliable. While 
the latter point is irrelevant for systems where the quantum ground state is close to 
the reference state (i.e. in our model in case of well-pronounced Neel order) it becomes 
relevant if the quantum ground state is far from the reference state (i.e., in our model 
when Neel order breaks down). Hence the results for higher spin quantum numbers 
must be taken with extra care. 
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Figure 4. Extrapolated sublattice magnetization M/s versus J' for various spin 
quantum numbers s using coupled cluster method (CCM), see text. 



We have calculated CCM results within the SUBn— n approximation for n = 2, 4, 6 
for s = 1/2,1,3/2,2. For spin 1/2 also results for n = 8 are available (see Refs. 
[TEll2"2| ). First we report the values for J' c for spin s — 1/2 (see also Refs. [HI 1221)- The 
extrapolation of the SUBn— n data for M with n = 2, 4, 6 as described above leads to 
J' c ~ 3.5. However, as discussed above the extrapolation of the order parameter tends 
to overestimate J' c (note that J' c obtained this way is even larger than the value found 
within MFA) and the extrapolation of the inflection point is favorable. We found as 
inflection points of the M S (J') curves J' inf {n) = 3.60 (SUB2-2), 3.33 (SUB4-4), 3.13 
(SUB6-6) leading to an extrapolated value of J' c = J' in Aoo) = 2.56. We mention that the 
consideration of SUB8-8 data leads to a slight modification of J' c to J' c = J' in f(oc) = 2.54, 
only. 

We now consider the case s — 1, where the results for the order parameter M are 
giving in Fig. |3J Clearly we see the weakening of the magnetic order by increasing 
J'. The extrapolation of the SUBn— n data for M with n = 2,4,6 leads to J' c ~ 11.7, 
i.e. we get the same tendency as for the variational MFA that J' c increases with s. 
Again the extrapolation of the order parameter leads to an overestimation of J' c . This 
overestimation is connected with the change in the sign of curvature of M(J') seen in 
Fig. El The favorable extrapolation of the inflection points leads to J' c = J' in j(oo) ~ 6.4, 
where the inflection points for the different levels of SUBn — n approximations are 
J' in = 2) = 3.93, J' (n = 4) = 6.04, J'(n = 6) = 6.36. 
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Table 1. CCM results for the ground state of the Heisenberg antiferromagnet on 
the square lattice with spin quantum number s — 3/2 and s = 2 using the SUBn— n 
approximation with n = {2,4,6}. Note that Np indicates the number of fundamental 
clusters at each level of approximation. For comparison we present the results of the 
second-order spin wave theory (SWT) [3*T]. 





S = 3/2 


S = 2 




N F 


E 9 /N 


M/S 


N F 


E 9 /N 


M/S 


SUB2-2 


1 


-4.943927 


0.936174 


1 


-8.593510 


0.950368 


SUB4-4 


15 


-4.976427 


0.910266 


15 


-8.633108 


0.93109 


SUB6-6 


461 


-4.982685 


0.89816 


461 


-8.640356 


0.922284 


SUBoo 




-4.9878 


0.8687 




-8.6461 


0.9011 


SWT 




-4.9862 


0.8692 




-8.6442 


0.9018 



Finally, we consider spin s = 3/2 and s = 2. The results for the extrapolated 
sublattice magnetization for spin values s = 1/2, 1, 3/2, 2 using SUBn— n approximation 
for n = 2, 4, 6 are shown in Fig.0J Evidently the sublattice magnetization Mj s increases 
with s demonstrating the decreasing influence of quantum fluctuations with growing spin 
quantum number. The critical value for s = 3/2 is obtained as J' c ~ 18.5 which is again 
too large in comparison to the MFA result. The extrapolation of the inflection points 
leads to J' c « 10.9. Note that we have calculated M using the CCM up to J' = 100 
for s = 2. However, we were unable to find a vanishing M (i.e., the critical value J' c 
obtained by extrapolation of the order parameter would be larger than 100.) Results 
for the point of inflection of M were similarly contradictory, and so the results for the 
position of the phase transition point predicted by the CCM for s = 2 are not included 
here. 

We conclude that the CCM SUBn— n approximation is inappropriate in order to 
describe the quantum phase transition correctly for higher spin values (namely, s > 3/2) 
at the levels of approximation currently available for present-day computers. However, 
we do observe that the tendency for critical value J' c to increase with growing s is 
observed using the CCM for s < 3/2, as expected. This problem of reliability might 
be resolved by going to higher orders of truncation index m, although we note that the 
computational problem is very difficult (e.g., with Np = 108033 for SUB8-8 for s = 3/2) 
and so this is not considered here. We note that LSUBn calculations do not place a 
restriction on the total number of spin flips used in the CCM correlation operators, 
although the fundamental clusters are restricted to remain within a locale defined by n. 
However, this again leads to an extremely large number of fundamental clusters even for 
low values of n and for higher spin quantum numbers, and so LSUBn is not considered 
here. Mean-field model states (e.g., based on the variational states in Sec. 3) might also 
provide enhanced results for the CCM. 
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As a byproduct we also present in Table [T] the results for the sublattice 
magnetization M/s for higher spin values for the pure square lattices (J 7 = 1), 
which are so far not calculated within CCM. We point out that for the pure square 
lattice the results for M are expected to be quite reliable, since the true quantum 
ground state is close to the reference state used as starting point. This is indeed 
confirmed by comparison with high precision second-order spin wave results [31] also 
presented in Table [0 We mention that due the reduced symmetry the number 
of fundamental configurations increases in case of J 1 ^ J. For SUB6 — 6 we find 
Np = 267, 1420, 1744, 1744 for s = 1/2, 1, 3/2, 2, respectively Note that N F for s = 3/2 
and s = 2 is equal only up to SUB6— 6 but differs for higher levels of approximation. 
For completeness we also give the sublattice magnetization for s = 1/2: M/s = 0.63 
(note that this value can be improved by considering also SUB8— 8 for the extrapolation, 
which yields M/s = 0.62 |2H]) and for s = 1: M = 0.81 (see also |3n|). 




Figure 5. Sublattcie magnetization M/s versus J' for spin quantum number s = 1 
using exact diagonalization of finite lattices of N — 8, 10, 16, see text. 



5. Exact diagonalization (ED) 

In addition to the MFA and the CCM we use the exact diagonalization Lanczos technique 
to calculate the order parameter for finite square lattices with periodical boundary 
conditions. The calculations are performed using the J.Schulenburg's spinpack |H2]- As 
usual for ED (see, e.g. Ref. 0) we calculate the square of the sublattice magnetization 
M 2 defined by M 2 = ([-^ J2f=i T i s i} 2 ) with the staggered factor r i( z A = +1, r i6B = — 1. 
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For the finite-size scaling of M 2 we use the standard three-parameter formula |3*3*1 I3"H |3] 
M 2 (N) = M 2 (oo) + dN- 1 / 2 + c 2 N~ l . The critical value J' c is that point where M 2 (oo) 
vanishes. Again we are faced with the problem, that the method becomes less reliable 
for larger quantum numbers s. While for s = 1/2 one can calculate the GS for the 
J- J' model up to N = 32 [TBI EI3 E2| sites one is restricted to lattices of up to N = 20 
for s — 1, up to N — 16 for s = 3/2 and up to iV = 10 for s = 3/2. Since for 
s = 2 we have only two lattices (N = 10 and N = 8) with the full lattice symmetry, 
we do not consider s = 2 within ED. To treat all three cases in a consistent way we 
consider only N = 8, 10, 16 for s = 1/2, 1,3/2. It is clear that the resulting finite-size 
extrapolation remains quite poor and allows only some qualitative conclusions. We 
present for illustration the results for the order parameter M for s = 1 in Fig. |S| The 
critical values obtained by finite size extrapolation of M are: J' c ~ 2.2 for s = 1/2, 
J' c ~ 5.5 for s = 1 and J' c ~ 10.1 for s = 3/2. These data confirm the tendency found 
by MFA and CCM that the increase of J' c with s is stronger than linear. 




Figure 6. The critical value J' c versus spin quantum number s obtained by different 
methods. MFA: variational mean- field approach (see Sec. |3J); CCM I: coupled cluster 
method (extrapolation of the order parameter, see Sec. @}; CCM II: coupled cluster 
method (extrapolation of the inflection point, see Sec. ED: exact diagonalization 
(see Sec. 
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6. Summary and Discussion 

We have investigated the ground-state magnetic order parameter for the square-lattice 
isotropic Heisenberg antiferromagnet with two kinds of nearest-neighbor exchange bonds 
(J- J' model) by using a variational mean-field approach (MFA), the coupled cluster 
method (CCM) and exact diagonalizations (ED). In particular, we have studied the 
influence of the spin quantum number s on the quantum critical point J' c . Our results 
for J' c are presented in FigEH and we note that a transition from a semi-classically Neel 
ordered phase to a magnetically disordered phase occurs at this point. Obviously, there is 
an increase of J' c with s signaling the diminishing of quantum effects. We have presented 
evidence that the critical value J' c increases with growing s according to J' c oc s(s + 1). 
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